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ESTIMATION OF SPATIAL CAPTURE DISTRIBUTIONS 
IN RESONANCE ABSORBERS 
by Donald Bogart 
Lewis Research Center 

SUMMARY 

The problem of precise calculation of spatial distributions of capture in resonance 
absorbers is crucial to the design of layered shields. Errors in spatial distribution of 
capture occur in multigroup neutron transport calculations because of the necessarily 
broad energy groups employed. The single average capture cross section in each group 
results in large underestimates of the capture rates near surfaces of resonance absorb- 
ers. Consequently, the spatial capture gamma ray generation and escape fraction are 
also in error. 

A method is presented for computing spatial resonance capture rates in thick layers. 
It employs group effective resonance integrals to precalculate group effective resonance 
cross sections that are universal functions of distance into the absorptive layer. The 
method is illustrated for captures in uranium-238 for the energy region 0. 5 eV to 
100 keV. 

The method is applied to a spherical reactor -shield configuration that contains al- 
ternate layers of depleted uranium and lithium hydride. Detailed comparison is made of 
the results of a discrete ordinates multigroup calculation with those of the present 
method. The comparison shows that the difference in spatial capture distribution of the 
S n broad treatment of resonance capture causes the capture gamma ray dose to be al- 
ways underestimated. For example, the difference in spatial capture distribution in a 
7-centimeter slab or uranium-238 causes the leakage dose to be a factor of 2 smaller 
than that obtained with the present method. The apparent generality of the present 
method suggests that it may be applied directly to the results of layered shield calcula- 
tions made by S n broad group methods. 

Application of the method to the experimental variation of epi-cadmium capture with 
depth from the surface of metallic uranium rods up to 5 centimeters in diameter as 
measured by Hellstrand provided spatial capture rates that agreed with the experiments 
very well. 


INTRODUCTION 


The problem of calculating the spatial distributions of capture in resonance ab- 
sorbers has not been given much attention because in the past only integral capture 
rates were of primary importance for reactor design. However, in the conceptual de- 
sign of shields for space power systems, the spatial distributions of capture gamma 
rays generated in metal layers play a crucial part in the shielding problem. 

Some attempts have been made to calculate spatial distributions of captures in 
resonance absorbers. For uranium rods the problem was approached by Richtmyer 
(ref. 1) as an illustration of the application of Monte Carlo methods. A spherical har- 
monics calculation of the radial distribution of neutron capture in a natural uranium 
rod in the NRX reactor was made by Kushneriuk (ref. 2) for the purpose of understand- 
ing the reactivity changes in fuel rods. 

The objective of the shielding calculation for a space reactor power source is to 
solve for the best arrangement of a heavy metal layered within a hydrogenous neutron 
attenuating material such that a given dose constraint is satisfied with a minimum shield 
weight. Metals such as tungsten and depleted uranium are considered to be advanta- 
geous even though resonance capture generates a considerable secondary gamma ray 
dose. 

Present shield design studies are made using multigroup neutron transport calcu- 
lations, but the question arises of how to handle the resonance region for thick layers 
of resonance absorbers. Because of digital computer memory storage limitations the 
calculations must be made with broad neutron energy groups each encompassing many 
resonances. In a recent study Lahti (ref. 3) evaluated multigroup and single broad group 
calculations of captures rates in a slab of depleted uranium. He investigated the 3 to 
11 eV neutron range that includes the large resonance at 6.67 eV, using the GAROL code 
(ref. 4) to evaluate group average cross sections. He found that the broad group failed 
to represent the multigroup spatial distribution, although the total capture rate was 
preserved. Therefore, the broad group calculation results in errors in the magnitude 
of the capture gamma rays that escape from the uranium slab. 

The object of this report is to present a method for the accurate estimation of spa- 
tial capture distributions in resonance absorbers. The method is illustrated for slabs 
of uranium-238 for six neutron energy groups covering the range from 0. 5 eV to 100 keV. 

The method is applied to a spherical reactor- shield configuration that contains al- 
ternate layers of depleted uranium and lithium hydride. Detailed comparison is made 
of spatial capture distributions from a discrete ordinates multigroup calculation with 
spatial capture distributions obtained by the present method. These comparisons illus- 
trate the effects of broad neutron grouping and the magnitude of resulting differences in 
capture gamma ray leakage due to incorrect calculation of spatial distribution of cap- 
tures. 
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The method is also applied to the experiments of Hellstrand (ref. 5) who measured 
the variation of epi-cadmium capture with depth from the surface of metallic uranium 
rods up to 5 centimeters in diameter. 


METHODS 

The procedure used here applies the results of integral resonance capture theory 
to estimate spatial distributions of captures within an absorber. Since the absorber 
will generally be surrounded by moderator regions with slowing down densities of dif- 
ferent magnitudes, the procedure will permit the estimation of the spatial distribution 
for a resonance absorber region in a nonuniform flux environment. 

Effective resonance integral. - The effective resonance integral is defined as 
that lethargy integrated absorption cross section per atom of absorber that will provide 
the correct capture rate C in an absorber of atom density N and volume V: 


C = NVI eff <p 


( 1 ) 


The flux (p is the scalar flux that exists in the moderating region prior to the in- 
sertion of the absorber. It is implicit in equation (1) that for nonabsorptive moderators, 
the flux varies inversely with neutron energy E. In order for equation (1) to be used 
quantitatively for a wide range of absorber size and geometry, the quantity I e ^£ must 
be calculated with precision. The many approaches to the problem of calculating I g ££ 
are summarized by Weinberg and Wigner (ref. 6), Dresner (ref. 7), and by Nordheim 
(ref. 8). It is the method of Nordheim, programmed for digital computer by Kuncir 
(ref. 9), and extended by Stevens (ref. 4) that obtains direct numerical solutions for the 
average fluxes in a resonance absorber surrounded by a moderator. The unresolved 
s-wave resonances are treated by the TUZ code of Kuncir. 

In the Nordheim method the average fluxes in the two adjacent media are obtained by 
solving the coupled integral equations for energy dependent fluxes that express the neu- 
tron balance in each region. The spatial transfers of neutrons between the adjacent 
media are effected by means of their respective collision escape probabilities and appli- 
cation of a reciprocity condition. A space-independent source is assumed in each re- 
gion for the calculation of these collision escape probabilities for spherical, cylindrical, 
and slab geometries. 

Remarkable success has been achieved in analytical confirmation of experimental 
values of I e ££ for many isotopes and absorbers of various geometry (see refs. 6 to 8). 
For a given isotope of various geometry, values of I eff have been correlated by a 
surface-to-mass parameter S/M, and it has become conventional to report effective 
resonance integrals as a function of this parameter. Dresner (ref. 7) has shown that 
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the average escape probabilities, that are required in the calculation of I e jj, are very 
nearly the same for infinite cylinders, slabs, and spheres when compared on the non- 
dimensional basis of the ratio of mean chord length l to potential scattering mean free 
paths. The value of l is given by 4V/S where V/S is the representative volume to 
surface ratio for the absorber. Therefore, the correlation of measured and calculated 
values of I e ^ for a variety of shapes using the surface-to-mass parameter S/M is a 
reasonable expectation. It is this correlation that is the basis for assuming that the 
present method that is evolved for slabs, may also be applied to other geometries. 

Effective resonance cross section . - The present method for estimating spatial dis- 
tribution of captures is derived for the case of a slab of thickness t that is infinite in 
the other directions. The slab is sufficiently thick so as to separate the moderating re- 
gions on both sides which in general support different slowing down densities and there- 
fore fluxes of different magnitude. However, for the slab in a symmetrical flux, we 
make use of the fact that half of the captures are due to the flux in the left region and 
half due to the flux in the right region. For the unsymmetrical case it is reasonable to 
expect the captures to divide in proportion to these surface region fluxes and to be inde- 
pendent of each other. Therefore, equation (1) may be written as the sum of the left and 
right capture rates per unit surface area: 


C = 


C L + C R 




( 2 ) 


where for the slab in a uniform flux, <p L and cp-^ are the equal scalar fluxes in the left 
and right moderating regions. 

Separating the left and right capture rates per unit surface area gives 


C 


L “ 





(3) 


The value of the product tI eff /2 as a function of slab thickness is proportional to 
the capture rate in the slab surrounded by a 1/E slowing down flux. The product 
approaches zero as t approaches zero and I e ^ assumes its maximum value of the 
dilute resonance integral. The product monotonically increases as t gets larger and 
I e ff assumes smaller and smaller values. 

A new function, which varies uniquely with distance x from the surface of the slab 
is defined in the next equation in which the left capture rates are differentiated with re- 
spect to x: 
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( 4 ) 


/tl 


dCj- 

dx 


eff ] 


= N<p, 


dt 


s N ^L a eff (x) 


Equation (4) defines a spatially dependent effective resonance cross section a g ^(x) 
that is a function of distance x from the left surface and expresses the change in slab 
capture rate with increment in slab thickness due to the scalar flux in the moderating 
region on the left. The variation of a e ff( x ) with distance x is obtained by numerical 
differentiation using calculated values of 1^. The variation of cr e ff( x ) with x varies 
from one-half the dilute resonance integral at the surface and falls rapidly with distance 
into the absorber. 

For the general case of a slab with different fluxes existing at its face, the total 
captures per unit surface area due to the specific fluxes in the left and right surface re- 
gions are then summed using the generalized values of a e ff( x ): 


«■/ 

•'0 


f 

» n 


N ^L a eff^ dx + / N<p R a eff (t - x)dx 


(5) 


EFFECTIVE RESONANCE INTEGRALS AND CROSS SECTIONS FOR URANIUM-238 

For general use, values of and have been calculated for uranium-238. 

The resonance region above 0. 5 eV has been divided in the resolved region to 2 keV into 
five groups with one additional group for the unresolved region extending to 100 keV. A 
number of groups are required to take care of moderators with a slowing down flux which 
varies in a manner other than 1/E. In general, a non- 1/E slowing down flux is present 
in a hydrogenous shielding material such as lithium hydride. It is implicit in any mul- 
tigroup calculation that a 1/E flux exists within each neutron energy group. 

The GAROL code (ref. 4) has been used to evaluate effective resonance integrals 
for the resolved region using the recommended resonance parameters (ref. 10) for 
uranium-238. Calculations were made for a wide range of slab thicknesses immersed 
in an infinite sea of hydrogenous moderator. For the unresolved region from 2 to 100 
keV the s-wave contributions were calculated using average s-wave resonance parame- 
ters determined from the resolved region and the code TUZ (ref. 9). Contributions to 
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TABLE I. - EFFECTIVE RESONANCE INTEGRALS FOR METALLIC URANIUM-238 


Slab thick- 
ness or 
rod radius, 
cm 

Geometry 

parameter, 

(S/M) 1 / 2 

Neutron energy group 

0. 5 - 10 
eV 

10 - 25 
eV 

25 - 100 
eV 

100 - 500 
eV 

0. 5 - 2 
keV 

2 - 100 
keV 

0. 5 eV - 
100 keV 




Effective resonance integral, I ef £ 

b 


1x10”® 

325.0 

133.65 

59. 84 

55.31 

17.40 

3.59 

2.78 

272. 57 

lxlO' 5 

102.8 

132.77 

59.44 

54. 95 

17.36 

3. 59 

2.78 

270. 89 

lxlO” 4 

32.50 

127.09 

56.87 

52. 56 

17.08 

3.59 

2.78 

259. 97 

5X10' 4 

14. 54 

111.83 

49. 98 

46.25 

16.26 

3.56 

2.77 

230.65 

.001 

10.28 

99.36 

44. 34 

41. 14 

15.48 

3.52 

2.76 

206.60 

.005 

4. 595 

57.29 

25.28 

24.47 

12.04 

3.32 

2.75 

125. 15 

.01 

3.250 

39.77 

17.39 

17.71 

9.82 

3. 13 

2.74 

90. 56 

.05 

1.454 

16. 10 

6.81 

8. 00 

4.92 

2.41 

2.65 

40.89 

. 1 

1.028 

11. 24 

4.66 

5.79 

3.58 

1.98 

2. 54 

29.79 

.5 

.4595 

5.25 

2.04 

2.87 

1.86 

1. 10 

2.22 

15. 34 

1.0 

. 3250 

3.86 

1. 44 

2. 15 

1.45 

.87 

2.06 

11.83 

2.0 

.2300 

2.85 

1.02 

1.61 

1. 15 

.71 

1.92 

9.26 

5.0 

. 1454 

1.88 

.63 

1. 07 

.90 

.56 

1.75 

6.79 

7.0 

. 1230 

1.60 

. 53 

.90 

.77 

! .51 

1.68 

5. 99 

10.0 

. 1028 

1.33 

.44 

.73 

.67 

.47 

1.61 

5.25 

15.0 

.0840 

1.06 

.36 

.55 

.57 

.41 

1.52 

4.47 

20.0 

.0727 

.89 

.32 

.46 

.50 

.38 

1.45 

mm 


the dilute resonance integral due to resonances other than s-wave were obtained by leth- 
argy integration of the recommended total capture cross section curve (ref. 10). 

Values of I e ^ are listed in table I and are shown in figure 1 as a function of the 
geometry parameter (S/M)^^. The dilute or unself -shielded integrals are approached 
for thin samples on the right and highly self-shielded slabs up to 20 centimeters thick 
are on the left. Such thick slabs are of interest to the shield designer and it should be 
noted that all resonance groups are important in the thick region. 

Numerical differentiation of each of these curves according to equation (4) provides 
the variation of the effective resonance cross section cr e ^(x) with the distance x from 
a surface of the slab. These results are listed in table II and are shown in figure 2 for 
each neutron group. Values of cr e ^(x) vary from one-half the dilute resonance integral 
at the surface to highly self- shielded values deep within a slab. Physically o" e ff(x) ex- 
presses the change in slab capture rate with increment in slab thickness due to the 1/E 
slowing down flux in the moderating region at one of its surfaces. The rapid reduction 
of Og ff (x) in the first few thousandths of a centimeter is due to the depletion of the avail- 
able flux at the resonance peaks; further reduction near the surface is caused by flux 
depletion in the wings of resonances. The more gradual reduction of <J eff (x). over the 
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TABLE IL - EFFECTIVE RESONANCE CROSS SECTIONS FOR 
METALLIC URANIUM-238 
ince Neutron energy group 


0.5-10 

eV 

10 - 25 
eV 

25 - 100 
eV 

100 - 500 
eV 

0.5-2 

keV 

2 - 100 
keV 

0. 5 eV - 
100 keV 


Effective resonance cross section, a 




1. 75 
1.72 
1.67 

1.59 1.37 

1.52 1.37 

1.45 1.36 

1.40 1.35 

1.23 1.34 

1.00 1.29 

.855 1.25 

.755 1.22 

.680 1. 19 

.483 1.11 

. 375 1.03 

. 337 .980 






Distance from surface, cm 

Figure 2. - Generalized group effective resonance cross sections for metallic uranium-238. 










bulk of the absorber is caused by flux attenuation of the low cross sections in the over- 
lap region between resonances. 

I 

Example. - As mentioned before, it is possible to estimate the distribution of spa- 
tial capture rate for thick slabs surrounded by moderators with different slowing down 
densities by independent consideration of group fluxes at the two surfaces of the slab. 
Such unsymmetrical spatial capture occurs in thick slabs of resonance absorbers con- 
sidered for layered shields. The present generalized slab procedure may be applied di- 
rectly to the results of broad group transport calculations and requires only the group 
fluxes at surfaces of absorber regions. 

As an example, the method is applied to a 4-centimeter-thick uranium slab for 
neutrons in the group 0. 5 to 10 eV. The left region flux <p L is maintained as unity 
while the right region flux <p R has values of (p R = 0 (corresponding to an outer shield 
layer of uranium), or (p R = <p L / 4 (approximating the flux falloff in a shield layer sur- 
rounded by moderator), or <p R = <p L (corresponding to a slab in a symmetrical flux). 
The distributions for these cases are shown in figure 3. The value of at the left 

surface is one-half the dilute resonance integral of 134 barns for 0. 5 to 10 eV neutrons. 
Also shown are the fractions of total capture up to several thicknesses. Relatively 



Figure 3. - Spatial distribution of captures per unit thickness in 4-centimeter-thick uranium slab for neutrons 
from 0.5 to 10 eV. 
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large capture fractions occur near the surfaces for slabs surrounded by moderator so 
that spatial capture gamma ray generation and escape from slabs are quite different 
from that of a transport calculation using a single group average capture cross section. 

Specific comparisons of doses due to gamma rays generated in slabs by various 
multigroup treatment for symmetric fluxes were illustrated by Lahti (ref. 3). Additional 
examples for slabs in flux gradients will be presented in the next section of this paper. 

A comparison of the present analytical method with experiments for rods is presented 
in the final section. 


APPLICATION TO LAYERED SHIELDS 

Applications of the discrete ordinates method to the calculation of layered shields 
for reactor sources are being studied at Oak Ridge (ref. 11) and at NASA (refs. 12 
and 13). At NASA (ref. 13), the weights of layered spherical shields for a fast spectrum 
nuclear reactor for space power have been calculated. The shield contains alternate 
layers of depleted uranium and natural lithium hydride and the thicknesses of the heavy 
layers were optimized (ref. 12) to meet a dose rate constraint of 2 millirem per hour 
at 20 meters from the shield. 

The one-dimensional discrete ordinates transport program ANISN (ref. 14) was 
used with Sjg quadrature and Pg cross sections to obtain neutron flux distributions 
for 26 energy groups throughout the spherical reactor- shield configuration. This par- 
ticular transport calculation is used to illustrate the differences in the spatial distribu- 
tions within the depleted uranium layers that result from representing the resonance 
energy region by the multigroup method and by the present method. For this particular 
problem, group average cross sections for the resonance region were also generated by 
the code GAROL and are part of the set of cross sections precalculated and tabulated for 
the convenience of shield calculations by Lahti and Westfall (ref. 15). 

The spherical geometry consisted of a reactor core, a pressure vessel, and a mo- 
lybdenum reflector that is surrounded by a seven layered lithium hydride shield contain- 
ing three layers of depleted uranium of thicknesses 6. 12, 4. 27, and 2. 90 centimeters. 
The spherical layers are sufficiently large in radius so that they may be considered to 
be equivalent to infinite slabs. 

The present method of estimating spatial capture distribution was applied to each 
layer of depleted uranium and for all neutron groups of the resonance region. It is noted 
that although the group GAROL average cross sections that represent the layered regions 
in the transport calculations can conserve total capture rates for the uranium layers, 
they do not produce the correct spatial capture rates. By a detailed comparison of the 
multibroad group transport results for this spherical shield with spatial capture results 
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generated by applying the present method, it is possible to calculate the differences in 
the capture gamma ray escape fraction for the two spatial distributions. These calcu- 
lations have been performed and the results for a few cases are illustrated. 

Spatial capture differences . - A detailed comparison of the spatial capture rates in 
the three uranium layers are shown in figure 4 for neutrons in the group from 25 to 
100 eV. The flux gradient across each uranium layer is primarily a function of its 
thickness and each thickness demonstrates a representative departure of the broad-group 
representation from a more accurate detailed spatial capture rate. In each case the 
broad one-group distribution shows essentially a single exponential fall- off of capture 
rate whereas the present method shows high surface region capture rates with rapid > 
falloff and relatively flat capture in the bulk of the layers. Although the total capture^ 
have been made identical by normalization of both distributions, it is the large differ- 
ences in capture rates at the outer surface regions of the layers that result in large dif- 
ferences in capture gamma ray escape. 

For each spatial distribution, a group dose factor DF^ to a detector at the right 
surface of each slab was computed for the uranium layers. Each uranium slab of thick- 
ness t was taken to be approximated by an infinite plane gamma-ray source of thickness 
dx and DF^ was computed as: 

f ^method® 1^ 

DF^t) = JL (6) 

f SWl-group 15 !^ 

J 0 


where S(x) are local values of captures per unit thickness, Ej(b) is the exponential in- 

/»oo 

tegral / (e~ a /a)da, and b is the product of the mass absorption coefficient at 


~2 MeV gamma-ray energy, uranium density, and distance x of source plane dx from 
the right surface of the uranium layers. These gamma-ray doses represent the uncol- 
lided components only. The dependence of the neutron capture gamma -ray spectrum on 
neutron energy may be important (ref. 16) in this calculation but has not been included 
in this example. 

Each group dose factor was in turn weighted by the fraction of total resonance cap- 
ture occurring in each group so as to arrive at a spectrally weighted dose factor for the 
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given slab layer thickness: 



( 7 ) 


The results of these calculations are presented in table III for the resonance region 
from 0. 5 eV to 100 keV which in general encompassed about 85 percent of the total cap- 
tures for the slab layers. The fraction of total captures and the dose factor are shown 
for six neutron groups in the resonance region from 0. 5 eV to 100 keV for which the 
present method was formulated. The results for the three thicknesses of depleted uran- 
ium show the spectrally weighted dose factor to increase with layer thickness. This re- 
sult is presented in figure 5 which shows the calculated dose factors to increase 


TABLE m. - RELATIVE DOSES DUE TO CAPTURE GAMMA RAYS IN 
SPHERICAL LAYERED SHIELDS OF DEPLETED URANIUM AND 
LITHIUM HYDRIDE FROM SPATIAL DISTRIBUTIONS OF 
CAPTURES CALCULATED BY PRESENT METHOD 
AND BY A MULTIGROUP DISCRETE 


ORDINATES METHOD 


Neutron energy 
group 

Thickness of depleted uranium layer, cm 

2.90 

4. 27 

6. 12 

Fraction 
of total 
captures 

Dose 

factor 

Fraction 
of total 
captures 

Dose 

factor 

Fraction 
of total 
captures 

Dose 

factor 

0. 5 - 10 eV 

0. 156 

1.32 

0. 143 

1.48 

0. 156 

1.81 

10 - 25 eV 

.079 

1.36 

.073 

1. 55 

.080 

1.85 

25 - 100 eV 

. 125 

1.33 

. 118 

1. 52 

. 135 

1.76 

100 to 500 eV 

. 132 

1.31 

. 138 

1.71 

. 138 

1.94 

0.5-2 keV 

. 237 

1.41 

.239 

1.61 

.237 

2.07 

2 - 100 keV 

. 105 

1 1.28 

. 115 

1.51 

. 114 

1.42 

0. 5 eV - 100 keV 

0.834 

1.35 

0> 826 

1.57 

0.860 

1.85 


14 




Figure 5. - Spectrum weighted dose factors for capture 
gamma rays in layered shields of depleted uranium and 
and lithium hydride. 


smoothly with layer thickness indicating a dose factor of 2. 0 for a uranium layer thick- 
ness of 7 centimeters. 

It is the large differences in calculated capture rates at the right surface regions 
that account for the large underestimate in capture gamma ray dose to a detector at the 
right surface of each uranium layer. The comparably large differences in calculated 
capture rates at the left surface regions do not affect the gamma dose at the right sur- 
face since these gamma rays are largely self-attenuated in traversing the uranium layer. 

Because of the apparent generality of the present method, it is suggested that the 
method may be applied directly to the results of discrete ordinates multigroup calcula- 
tions in order to estimate more accurately the spatial distributions of captures in thick 
slabs. 


APPLICATION TO EXPERIMENTS WITH RODS 

Some excellent experiments have been performed by Hellstrand (ref. 5) that meas- 
ured the spatial distributions of resonance captures in rods of uranium. Because of the 
rapid depletion of resonance flux near the surface, it may be expected that the universal 
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functions ^ e ff( x ) for slabs, shown in figure 2, will also apply for rods since the capture 
rate falls off rapidly at depths that are very much smaller than the radii of curvature of 
the surfaces. 

In the Hellstrand experiments, the variations of epi-cadmium capture with depth in 
metallic uranium rods up to 5 centimeters in diameter were measured. The spatial dis- 
tributions in rods were determined as part of a series of measurements of effective 
resonance integrals for uranium metal and oxide in various geometries. Resonance cap- 
ture was measured by the activity of neptunium-239 induced in samples and the surface 
resonance flux was monitored by activation of thin gold foils. The results were evalu- 
ated in barns by using irradiated thin uranium samples as comparative standards using 
a dilute resonance integral of 280 barns for infinitely thin samples. 

For these uranium rods, samples of metal at various depths were obtained by turn- 
ing central sections on a lathe and collecting and counting the cuttings. The uncertainty 
in depth location of layers of measured rod activity near the surface due to machining 
was claimed to be ±0. 00012 centimeter. It should be noted that in this arrangement, 
samples of activated uranium collected near the surface of rods are exposed effectively 
only to the flux near the surface. On the other hand, induced activities are evaluated 
relative to infinitely thin foils of uranium which are activated by the unshielded flux on 
both surfaces. Therefore it should be expected that the activity induced near the surface 
of samples will approach a value of 140 barns or one-half of the uranium dilute resonance 
integral. Hellstrand designated these relative spatial activities as effective resonance 
cross sections which shall be seen to correspond to the values of capture rate per atom 
evaluated by the present method. 

All cadmium covered rods were irradiated in the central channel of the Swedish 
heavy water reactor. In this connection it was pointed out by Hellstrand that all of his 
results do not include capture in the energy region from fission energy down to about 
100 keV. This is so because the reactor location where samples were irradiated was 
removed from the reactor fuel and only a small part of the incident flux spectra extended 
above 100 keV. Therefore the calculated effective resonance integrals for the range 
0. 50 eV to 100 keV shown in table I apply directly to these experiments except for a con- 
vention observed by Hellstrand. This convention was to report all of his results as epi- 
cadmium non-l/v effective integrals. Hellstrand assumed the 1/v part of the integral 
to be constant and equal to 1. 1 barns corresponding to a cadmium cutoff of about 0. 6 eV 
and a thermal absorption cross section of 2. 75 barns. In order to compare the present 
integral and spatial calculations with Hellstrand, a comparable 1/v part of the integral 
was subtracted from the GAROL results. 

A comparison of the calculated and measured epi-cadmium non-l/v effective reso- 
nance integrals for uranium rods is given in table IV. Over the range measured the 
calculated GAROL results are in excellent agreement with Hellstrand' s empirical fit to 
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TABLE IV. - COMPARISON OF CALCULATED 
AND MEASURED EPI-CADMIUM NON-l/V 
EFFECTIVE RESONANCE INTEGRALS 
FOR URANIUM IN THE RANGE 


0. 50 eV - 100 keV 


Slab thick- 
ness or rod 
radius, 
cm 

Geometry 

parameter, 

(S/M) 1 / 2 

Effective resonance 
integral, I eff , b 

GAROL 

Hellstrand 

(a) 

0.01 

3.250 

89.4 

83. 1 

.05 

1.454 

39.7 

38.7 

. 1 

1.028 

28.6 

28.2 

.5 

.4595 

14.2 

14.2 

1.0 

.3250 

10.8 

10.8 

2.0 

.2300 

8.3 

8.5 

5.0 

. 1454 

5.8 

6.4 

7.0 

. 1230 

5. 1 

5.8 

10.0 

. 1028 

4.4 

5.4 

15.0 

.0840 

3.6 

4.9 

20.0 

.0727 

3.2 

4.6 


a From Hellstrand's fit to experiment: 
I eff = 2.81 + 24.7(S/M) 1/2 ; 

0. 2 < (S/M) 1 / 2 < 0. 8. 


his experiments; however for larger and smaller rods, the GAROL results deviate from 
this limited empiricism. 

The question now arises of how to apply the values of a g ^(x) generated for the case 
of slabs to the case of infinitely long cylindrical rods of radius R placed in a moderator 
supporting a uniform slowing down flux. Although the variations of cr ef j(x) with x are 
derived for slabs we assume that they also apply to a good approximation for rods of 
resonance absorber. This is based on the surface-to-mass correlation of for rods 
and slabs for which the same values of I g ^ apply for equal values of surface-to-mass 
parameter (S/M) which generalizes results for many absorber geometries (see Dresner 
(ref. 7)). 

Atoms near the surface of rods will have values of Og^x) of one-half the dilute 
resonance integral while atoms near the center of the rods will have values of 
corresponding to the value of x equal to R; intermediate values of <7 e ff( x ) are assumed 
to correspond to the position along the diameter of the rod. With this assumption, the 
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Figure 6. - Comparison of spatial distribution of captures calculated by present method with experimental data of Hellstrand for metallic rods of 
uranium-238. 
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Figure 7. - Comparison of experimental and calculated spatial distribution of captures near surface of metallic 
rods of uranium-238. 


spatial variation across the entire rod diameter is evaluated in the same manner as for 
slabs in which the distance from the rod surface x varies from 0 to 2R. 

The total captures per unit length of rod is given as 



The present method for estimating spatial distributions of captures has been applied 
to the rods measured by Hellstrand using the GAROL epi-cadmium non-l/v effective 
resonance integrals from 0. 50 eV to 100 keV to provide generalized values of cr^. 
These values of cr e ^ are then used to calculate spatial capture distributions per atom 
that are directly comparable to the Hellstrand measurements. These comparisons are 
presented in figures 6 and 7. The calculations are absolute with no normalization and 
the agreement is particularly good over the bulk of the volume for all rods. In the small 
volume region near the surface the present method provides results in remarkable 
agreement with experiment. All the rods show the same rapid falloff of capture 
rate per atom from its value of one-half the dilute integral of 136 barns at the surface. 
Therefore, it may be concluded that the present method for estimating spatial distribu- 
tions of captures in uranium slabs predicts experimental spatial distributions of cap- 
tures in rods very well. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, January 14, 1970, 

124-09. 
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